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We study the two-qubit controlled-not gate operating on qubits encoded in the spin state of a 
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two quantum dots encoding a single qubit, while tunneling between the quantum dots that belong 
to different qubits is forbidden. Therefore, the two qubits interact exclusively through the direct 
Coulomb repulsion of the electrons. We find that entangling two-qubit gates can be performed by 
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quantum dots within the qubits. The entangling interaction can be controlled by tuning the bias 
through the resonance between the singly-occupied and doubly-occupied singlet ground states of a 
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I. INTRODUCTION 

The spin- 1/2 of a single electron trapped in a quantum 
dot (QD) is a promising candidate for a carrier of quan- 
tum information in a quantum computer Q . To perform 
a quantum computation we need to have all the unitary 
operations from some universal set of quantum gates at 
our disposal 2] . One such universal set consists of all the 
single qubit quantum gates and a two-qubit controlled- 
not (CNOT) quantum gate. Quantum computation over 
the single-spin qubits with the logical states correspond- 
ing to the spin orientations ||) and ||) can in principle 
be achieved using an external magnetic field or with g- 
factor engineering for the single qubit operations, and 
with the time-dependent isotropic exchange interaction 
H ex {t) = J(t)Si • S2 for manipulating a pair of qubits 
encoded into spins Si and S2 

Control of electron spins in quantum dots is in the 
focus of many intense experimental investigations. Ma- 
nipulation of pairs of electron spins using the tunable 
isotropic exchange interaction has already been demon- 
strated in several experiments 0, 0, 01 ■ Such control was 
used in a study of the QD spin decoherence due to the hy- 
perfine coupling to the surrounding nuclear spins, where 
the splitting between the singlet states with the total spin 
S = and the triplet states with S — 1 was used to turn 
on and off the singlet-triplet mixing caused by the hyper- 
fine interaction. An important result of these studies is 
that the coherence time of an electron spin in a quantum 
dot is very long if the decoherence due to the interaction 
with the nuclear spins can be suppressed. The spin co- 
herence times can be improved by the manipulation of 
nuclear spins 0, 0, @ j m principle allowing for elaborate 
sequences of operations to be performed. Single spin con- 
trol is based on the local manipulation of the magnetic 
field or g- factor Q] , or on ESR methods 0, 0] and has 
only recently been demonstrated experimentally |Tl) . 

The difficulty of single-spin control has inspired a num- 
ber of proposals for quantum computation based on the 
encoding of qubits into more than one spin. These encod- 
ing schemes reduce the requirement on the control over 



electron spins, but have the drawback of introducing so- 
called leakage errors in which the state of encoded qubit 
"leaks" out of the set of computational states. Standard 
error-correction procedures can be modified to prevent 
this kind of errors • A universal set of quantum gates 
operating on qubits encoded into states of three quantum 
dot spins with equal total spin quantum numbers can be 
implemented through control of the isotropic exchange 
coupling alone H ex [H Q Control over interac- 

tions that are symmetric only with respect to rotations 
about a fixed axis in spin space allows for the construc- 
tion of a universal set of quantum gates that operate over 
qubits encoded into pair of spins. One such encoding is 
into the orthogonal states |t|) and |||) of two spins-1/2. 
A universal set of quantum gates over such qubits can in 
principle be performed by the control over H ex , with the 
anisotropy provided by an external static homogeneous 
magnetic field and a site-dependent g-factor [0| [jji ■ 

We consider a variant of the two-spin encoding where 
the logical zero |0l) and the logical one quantum 
states are the singlet and the triplet with zero projection 
of the total spin to the symmetry axis z (S z — 0), e.g., 
for lateral QDs, the z axis is the normal to the plane of 
the heterostructure, 
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These qubits can be manipulated by an axially sym- 
metric interaction to produce a universal set of quan- 
tum gates. The interaction with an external magnetic 
field and the isotropic exchange [l6l Il7| , or the interac- 
tion with an external magnetic field and an anisotropic 
spin-orbit coupling [Tsj . or the spin-orbit coupling alone 
[19| , were all proposed as a way of producing a universal 
set of quantum gates operating on singlet-triplet two- 
spin qubit, Eq. JIJ. Recently, it was suggested that an 
architecture based on singlet-triplet qubits individually 
addressed using isotropic exchange interaction and inho- 
mogeneous magnetic field and coupled through Coulomb 
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FIG. 1: Double-double quantum dot (DDQD) setup. The 
four single-electron quantum dots are aligned along a fixed 
direction. The spins of the electrons on two quantum dots, 
inner (I) and outer (O), separated by a distance 2oi encode 
a qubit. Two such double quantum dot (DQD) qubits, left 
(L) and right (R), at the distance 2<i2 are separated by an 
impenetrable barrier. The tunneling matrix element t within 
the double quantum dots (DQD) carrying the qubits, and the 
bias e of the inner dots with respect to the outer are equal 
on both DQDs, and can be electrically tuned. The Coulomb 
interaction between the DQD is represented by the capacitor 
C. 

interaction of the electrons is scalable and in principle 
realizable |2fJ. 

In this paper, we study a particular realization of en- 
tangling two-qubit gates between singlet-triplet qubits, 
Eq. Pjl. where each qubit is represented by a pair of 
tunnel-coupled single-electron quantum dots, as pro- 
posed in |20| . In this realization, the double quantum 
dots are separated by a barrier which is impenetrable 
for the electrons, so that the qubits are coupled exclu- 
sively through the Coulomb repulsion of electrons, while 
the exchange terms between electrons on different double 
quantum dots vanish. The setup of this double-double 
quantum dot (DDQD) is illustrated in Fig. ^ 

The Coulomb interaction is spin-independent, lead- 
ing to an isotropic interaction JSi ■ S2 between tunnel- 
coupled spins Si and S2. The anisotropic correction to 
this interaction is dominated by the spin-orbit coupling 
induced term J/3 • (Si x S 2 ) + 0{\(3\ 2 ). The relative 
strength of the anisotropic interaction in the quantum 
dot systems in GaAs is estimated to be |/3| ~ 0.1 — 0.01 
[2TI I23 • The influence of the anisotropic corrections can 
be reduced in specific implementations of the quantum 
gates [23l . In our study of a two-qubit gate operation, 
we will only consider the case of isotropic interaction and 
neglect the weak anisotropy. In this case, transitions be- 
tween spin-singlet and spin-triplet states on a DQD are 
forbidden. Due to this spin symmetry, the four-electron 
Hamiltonian is block-diagonal, 

H =diag(H ss ,HsT,H TS ,H T T). (2) 

The non-zero blocks flab, where a, 6 — S,T, act on the 
states in which electron pairs on each DQD are either in 
the singlet (S) or in a triplet (T) state of the total spin 
5 = or 5 = 1. 



Our main results are the effective low-energy spin in- 
teraction and a scheme to perform a two-qubit CNOT 
gate in an electrically controlled DDQD system. The ef- 
fective low-energy spin interaction in this setup has the 
form 

H = J (S LI ■ S LO + S RI ■ S RO ) + E e \SS)(SS\. (3) 

Two pairs of spins, Sli and Slo on the left (L) qubit and 
Sri and Sro on the right (R) qubit (see Fig. inter- 
act via the isotropic exchange interaction of strength J, 
and the entangling interaction of strength E e that shifts 
the energy of the singlet-singlet state. We show how the 
entangling two-qubit quantum gates for universal quan- 
tum computation can be performed through the electrical 
control of E e . 

The triplet states with S z — 0, ±1 are degenerate in 
the absence of a magnetic field. A uniform magnetic field 
B, pointing along the z axis normal to the plane of QDs 
causes a Zeeman splitting g/i^B • S between the S z — 
states and the states with S = 1, S z = ±1. Our results 
apply both to the isotropic (B = 0) and anisotropic, but 
axially symmetric (B 7^ 0) case, if we take the S z = 
state to represent the qubit | Xjc) state. 

A two-qubit quantum gate can in principle be per- 
formed by adiabatically varying the tunneling matrix el- 
ement t and the bias e within the DQD. In practice, it is 
much simpler to change the bias e while t remains fixed 
[2?| . The control parameters e and t have to vary slowly 
on the time scale set by the energy splitting between the 
states of a given spin configuration. During the gate ap- 
plication, the orbital components of the S and T states 
are different due to the Pauli principle that forbids the 
electrons in a spin triplet to share their orbital state, see 
Fig. [5] As opposed to t and e that are determined by 
gate voltages and can be changed more or less at will, 
the Coulomb interaction is set by the geometry of the 
system and therefore fixed. We show how the control of 
the parameters t and e, or even e alone, can nevertheless 
be used to implement entangling two-qubit gates on en- 
coded singlet-triplet qubits through its influence on the 
Coulomb terms. 

When an adiabatic gate is applied, the lowest energy 
state in each block H a b, of energy Eat,, where a,b = S, T, 
see Eq. @ acquires a phase <j) a b = f*f E a b(t')dt' /K. The 
energy E a \y becomes time-dependent through the time- 
dependence of the parameters t and e in the interval ti < 
t' < tf. The resulting interaction is described by an 
effective 4-dimensional two-qubit Hamiltonian acting in 
the space spanned by the lowest-energy states \SS), \ST), 
\TS), and \TT) in the corresponding blocks H a b, and has 
the form of Eq. 10 . 

In the regime of strong bias, |e — U\ 3> t, where U is 
the on-site Coulomb repulsion, we investigate the DDQD 
system using perturbation theory. For the case of arbi- 
trary bias e, we numerically diagonalize the Hamiltonian 
Eq. J2J. We show that the two-qubit quantum gate can 
be operated by tuning the bias e so that the amplitude of 
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FIG. 2: Two-qubit quantum gate, a) When the inner quan- 
tum dots of the two double quantum dots system are strongly 
biased (e > U + t) the ground state is the doubly occu- 
pied inner dot. Due to the Pauli principle, only the spin 
singlets (5) can tunnel into the doubly occupied states on 
their DQDs. As the bias e is reduced, the states again be- 
come degenerate, b) A quantum gate is performed by send- 
ing a bias pulse e(t'). Each qubit state \ab) acquires a phase 
<t>ab = f_ E a b(t')dt' /h, where E a b(t') is the ground state en- 
ergy of the Hamiltonian at time t' reduced to the appropriate 
spin subspace, resulting in a two-qubit quantum gate. 

the doubly occupied state in the lowest energy spin sin- 
glet becomes appreciable. In this "on" state with a large 
double occupancy amplitude, entanglement is generated 
between the two-spin qubits. The entanglement genera- 
tion is suppressed in the "off" regime with weak bias and 
tunneling. Therefore, the generation of entanglement be- 
tween the two two-spin qubits encoded into DDQD can 
be efficiently controlled using the bias e alone. Together 
with the single-qubit operation this control is sufficient 
for universal quantum computing. 

This paper is organized as follows. In Sec. II, wc intro- 
duce our model of the DDQD system, followed by the dis- 
cussion of the control through voltage pulses. In Sec. Ill 
we focus on the case of the strongly biased (|e — U\ t) 
DDQD system and calculate the interaction between the 
qubits. The constraint of the strong bias is lifted in 
Sec. IV, where we numerically find the interaction be- 
tween the qubits, valid at an arbitrary bias e. In Sec. V, 
we outline the construction of a CNOT gate based on 
the resources for the control over a pair of qubits de- 
duced from the results of Sees. Ill and IV. Our results 
are summarized in Sec. VI. The technical details of the 
calculation are collected in the Appendix A. 



II. MODEL 

For the purpose of finding the effective low-energy spin 
Hamiltonian, the excited orbital states of single quan- 
tum dots can be neglected, leading to the Hund-Mullikcn 
(HM) approximation with one orbital per dot 0, |25| . 
In the HM approximation, the state space of the two- 
electron system in a double quantum dot (DQD) encod- 
ing the left (q = L) or the right (q — R) qubit is spanned 
by three singlet basis states, |5), \Di) and \Do) and one 
triplet basis state |T ) 

\~S) = ^ qI ^ q0l -c\ n cl m )n, (4) 

\ D o)=c\ o ^ qO[ \0), (6) 

l r o) = ^(4/T c Joi+ c Wor)|0)» ( 7 ) 

where Ck is the annihilation operator for an electron in 
the state k = {qk,Pk,Sk) on the qubit qu = L,R, with 
position pk = 1,0, where I stands for inner and O for 
outer quantum dot within a qubit, and spin =T-I- 
The vacuum |0) is the state of empty QDs. 

In the standard notation the singlet states of a DQD 
are denoted by \(n, m)S), where n is the number of elec- 
trons on the left QD and m is the number of electrons 
on the right QD. Our singly occupied singlet is then ex- 
pressed as |5) = |(1,1)5). The doubly occupied singlet 
states on the left, q = L, DQD are \Di) = |(0, 2)5), and 
\D ) = 1(2,0)5). On the right, q = R, DQD the defini- 
tions are reversed, \Di) = |(2,0)5), and \D ) = 1(0,2)5). 

The orbital states annihilated by c& approximate the 
ground states of the single-particle Hamiltonian 

^=E^(p*-~ A H) 2 +nr i ), (8) 

i 

describing an electron in the magnetic field B — V x 
A and confined to the system of quantum dots by the 
electrostatic potential V. The quantum dots form in the 
minima of this potential, which is locally harmonic with 
the frequency u)q. The ground states of Hi localized in 
these wells are the translated Fock-Darwin states jy] ■ 
The HM Hamiltonian is of the generic form 

H =t Yl ( 5 <lk,qi 5 Sk,si c l k Is k c 1kOs k +h.c.J - 
k.l 

1 , ( 9 ) 

e X! 4 c fc + 2 E ( kl \ v c\mn)clclc n c m . 

The intra-DQD tunneling term oc t preserves the electron 
spin. The bias e of the inner (p^ = I) QDs with respect to 
outer (j>k = O) is taken to be symmetric, i.e., the energy 
of both inner dots is lowered by the same amount. The 
two-body Coulomb interaction is denoted by Vc- Near 
the center of the quantum dot, the electrostatic potential 



is approximately harmonic and we assume that the wave 
functions of the electrons annihilated by the operators 
Ck are well approximated by the orthogonalized Fock- 
Darwin ground states. 

The impenetrable barrier that separates the DQDs im- 
poses the conservation of the number of L(R) electrons, 

= Ep=/,0; S =T,i" i («)P ;s ' Where ^PS = C \ps c qps- 

The fiLtm conserving terms, proportional to the interac- 
tion matrix elements (kl\Vc\mn) in Eq. JJJJ, where the in- 
dices k, I, to, n denote the single QD ground states, can be 
divided into intra-DQD terms where qf. = qi = q m = Qn 
and inter-DQD terms that satisfy qk ^ qi and q m 7^ q n . 
All the other terms, e.g., the ones that annihilate two 
electrons on the left (L) DQD and create two on the 
right (R) DQD violate the conservation of the electron 
numbers and therefore vanish. 



A. Interaction within a double quantum dot 

The terms for the interaction within a DQD in Eq. @ 
were discussed in 0. They renormalize the one-body 
tunneling matrix element t — > tjj = t+(S\Vc\Dj( )) /v2, 
introduce the on-site repulsion U = (-Dj(o) |Vc|Ar(o)) of 
two electrons on the same QD, and cause transitions be- 
tween the two doubly occupied DQD states with the ma- 
trix element X — (-Dj(o) | Vc?| ^O(i) )• Also, the Coulomb 
interaction on a DQD contributes V+ = (S\Vc\S) to 
the electrostatic energy of the symmetric and V— — 
(Tq\Vc\Tq) to the antisymmetric singly occupied orbitals 
of two electrons in a DQD [9j, giving rise to a direct ex- 
change interaction between spins. As a result, the elec- 
trons on a DQD are described by an extended Hubbard 
model with the isotropic exchange interaction |2( 



J = V- - V+ 



Uh 
2 



(10) 



where Uh = U—V+ + X is the effective on-site repulsion. 



B. Interaction between the double quantum dots 

The Coulomb interaction between the DQDs produces 
three new classes of direct terms in the Hamiltonian, 
while the exchange terms between the DQD vanish due 
to the impenetrable barrier. 

In the first class are the terms proportional to the num- 



ber operators h„ 



/, describing the electrostatic re- 



pulsion of the electrons in states qps and qp's', where 
L = R and R = L. For a pair of identical DQDs, there 
are three such terms: the interaction of a pair of electrons 
on the inner QDs, Un = (qIs,qIs'\Vc\qIs,qIs'}, the in- 
teraction of an electron on the inner QD of one DQD and 
an electron in the outer QD of the other DQD, Um = 
(qIs,qOs\Vc\qIs,qOs), and the interaction of electrons 
on the outer QDs, Up — (qO s , qO s\Vc\qO s , qO s) , 
Fig.EJi- 
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FIG. 3: Effects of the direct Coulomb interaction between 
double quantum dots (DQDs). All the exchange terms be- 
tween the DQDs vanish due to the impenetrable barrier, a) 
The Coulomb repulsion between the electrons on different 
double quantum dots contributes to the energy of the system. 
In the case of identical DQDs separated by an impenetra- 
ble barrier, there are three such contributions, coming from 
the electrons in orbitals that are near (Un), at a medium 
distance (Um) or far apart (Uf)- b) The tunneling matrix 
elements within a DQD are renormalized by Tj or To, due to 
the interaction with an electron on the inner or the outer dot 
of the other DQD. c) The interaction enables the correlated 
hopping processes in which electrons simultaneously tunnel in 
both DQDs. In one such process the electrons tunnel to the 
same side (either left or right) with the matrix element Xs- In 
the other correlated hopping process electrons simultaneously 
tunnel into the inner or outer quantum dots of their double 
quantum dots with the matrix element Xo- 



In the second class are the terms proportional to 

™<ips c qp'8' c <iP's' ' wnere 1 = and = 1. These terms 
describe the spin-independent correction to the tunnel- 
ing matrix element in the q qubit due to the interac- 
tion with an electron in the state qps. The two pa- 
rameters that determine the tunneling corrections are 
T p i = (qps, qp' s'\Vc\qps, qp' s') , and are due to the in- 
teraction with an electron in the p' = 1,0 orbital in the 
other DQD, Fig. Eb- 

The terms in the third class are proportional to 
c qps c qps c qp' s ' c qp's' i an d describe the processes in which 
electrons in both DQD tunnel simultaneously, Fig. Et- 
The two independent matrix elements for these pro- 
cesses are Xs = (qps,qps'\Vc\qps,qps') describing the 
tunneling from the inner to the outer orbital in one 
DQD and from the outer to the inner in the other, and 
Ad = (qps,qps'\Vc\qps,qps') describing the simultane- 
ous tunneling into inner or outer orbitals in both DQDs. 
For the system in zero magnetic field these two matrix 
elements are equal, A5 = Xr>. 



C. Control of the interaction 

In order to describe the influence of the intra-DQD tun- 
neling t and the bias e on the spectrum of the DDQD, 
we have to model the dependence of the Hamiltonian on 
these external parameters. In an experiment, both t and 
e are controlled by applying voltages to the electrodes 
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that define the quantum dots. The exact form of the 
voltage-dependent DDQD binding potential was studied 
using the Schrodinger-Poisson equation j2fj , but here we 
do not attempt to calculate the dependence of the Hamil- 
tonian Eq. J5J on e and t from first principles. 

Instead, we adopt a quartic double-well model for the 
potential of a DQD centered at (±02, 0) of the form |j| 



V(x,y) 



1 



((x T a 2 y 



y 



(li) 



where m is the electron effective mass, 2a\ is the distance 
between the approximately harmonic wells in a DQD, 
and 2a2 is the distance between the DQD double-well 
minima. In the limit of well separated dots, 3> aB, 
where ae is the QD Bohr radius given by eig = h/muio, 
and near the local minima of the quartic potential well 
at (±a2±ai, 0), the potential is approximately harmonic 
with the frequency luq. The Fock-Darwin ground state 
wave functions in this harmonic potential centered at 
{x c , 0) and in the magnetic field B normal to the plane of 
the dots, described in the symmetric gauge by the vector 
potential A = B(— y, x, 0)/2, are 



<t> Xc {x,y) = !I^ e -M(^^ 2 +v 2 )/^+"^^v/ r \ (12) 



where ujl — \J eB/2mc is the electron Larmor frequency 
and uj = y/ u>q + u)\ is the resulting confinement fre- 
quency with both electrostatic and magnetic contribu- 
tions. We will use the magnetic compression factor 
b = lu/ujq to measure the strength of the magnetic field, 
consistently with the notation in 

The translated single-electron Fock-Darwin states 
4>±a 2 ±a 1 {x,y) define the state space of the variational 
HM approximation for a DDQD. The tunneling matrix 
element between the Fock-Darwin ground states in the 
local minima of the potential Eq. is our control pa- 
rameter t la!, 



t = (^±O2+Ol|-ffl|0±O 2 -Ol) 



s 



8 1 + S* 2 



(13) 

where S = {(j}±a 2 + ai \<f>± a2 - ai ) = exp(-df(2b- 1/6)), is 
the overlap between the Fock-Darwin ground states in a 
DQD. 

As t is changed by external voltages, we assume that 
the overlap S between the oscillator states remains con- 
sistent with the relation Eq. (|13fl which is valid for the 
double-well potential V. All the Coulomb matrix ele- 
ments can be expressed in terms of S so that after solving 
equation l|13|) for the overlap they become functions of t, 
see Appendix A. The bias e is modeled as an energy shift 
of the orbitals, so that the inner pk — I orbitals have 
their energy reduced by e. 

The two-qubit gates are applied by time-dependent 
tuning of the tunneling matrix element t and/or the bias 
e in the DQDs using voltage pulses. In a typical experi- 
ment, the control of the QD energies through e is much 
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FIG. 4: Illustration of the double quantum dot energy levels 
as a function of the bias e. The energy of the singlet state 
with doubly occupied outer quantum dot, \Do), is indepen- 
dent of the bias. The energies of the singly occupied singlet, 
\S), and the singly occupied triplet, |To), state are lowered 
with the increasing bias as they have a contribution — 2e from 
the biased inner quantum dots. The energy of the singlet 
with doubly occupied inner quantum dots, \Di), is lowered 
with the increasing bias faster than the energy of \S) and \To) 
state, due to the bias contribution of — 4e. When the tunnel- 
ing t is zero, the lowest energy levels cross at the bias U±, 
leading to a drastic change of the effective spin interaction. 
For nonzero tunneling, the levels anticross, but the effective 
spin interaction still changes significantly when we tune the 
system from one side of the anticrossing to the other. 



easier to achieve than the control over tunneling matrix 
element t [2?! ■ The reason behind this is that the energy 
bias is linear in applied voltage, while the tunneling is 
typically exponential. 

The structure of the energy levels is particularly sim- 
ple in the limit of zero tunneling t — 0. In this limit, the 
eigenstates are the Hund-Mulliken basis states, Eqs. 
0. Their energies are determined by the bias e, the 
external magnetic field B, and the direct Coulomb in- 
teraction that is set by the device geometry. A drastic 
change in the structure of the DDQD spectrum as a func- 
tion of bias e appears at the crossings of the lowest energy 
singlet states within a DQD. Each of the singlet states 
\S), \Dj), and \Do) is lowest in energy for some values 
of the bias e, Fig. 0] A crossing occurs when either the 
positive bias overcomes the effective on-site repulsion U, 
making the state with both electrons in an inner dot \Dj) 
the lowest in energy, or the negative bias makes \Dq) the 
lowest in energy, see Fig. [5] We use the effective on-site 
repulsion U to emphasize the fact that it includes not 
only the repulsion of two electron in the same dot, de- 
noted by U, but also the energy of the interaction with 
the electrons on the other DQD. We will also use two 
special values of the effective on-site repulsion, U±. Due 
to the dependence of the effective on-site repulsion on the 
state of the other DQD, the lowest energy singlet-singlet 
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FIG. 5: Bias dependence of the double-double quantum dot 
(DDQD) ground state, (a) When the bias e of the inner quan- 
tum dots with respect to the outer ones is weaker than the 
effective on-site Coulomb repulsion U, the charge configura- 
tions of the lowest energy singlet and triplet states consists of 
singly occupied orbitals. (b) When e > U, the lowest energy 
singlet has a doubly occupied inner quantum dot, while the 
orbital state of the lowest energy triplet remains unchanged. 



DDQD state can consist of different singlets on the two 
dots, as in \S, Dj) and \Dj, S). In the strong bias regions, 
the lowest energy singlets are doubly occupied states. For 
e — U+ 3> t the lowest energy singlet is \DjDi), and for 
U- — e ^> t, the lowest energy singlet is \DoDo)- The 
second doubly occupied singlet state is separated by an 
energy gap « \2e\ from the lowest energy state. 



III. STRONG BIAS 

To develop an intuitive picture of the operation of an 
entangling two-qubit gate and the mechanisms for its 
control, we consider the simple case of strong bias. We 
show how the switching between the strong bias regime 
(e — U+ t), and the weak bias regime in which the dom- 
inant interaction is the on-site repulsion provides us with 
control over the entangling interaction E e . The boundary 
of the strong bias regime considered here is set by U+ = 
(3U N - 2U M -U F - 2V- + 2U)/2. A similar strong bias 
regime with the lowest energy singlet \Do, Do) exists for 
[/_-e> t, where U- = (3U F -2U M -U N -2V-+2U)/2, 
but we do not consider it here in detail. In both of these 
regimes, a wide energy gap « 2\e\ to the second doubly 
occupied state allows us to neglect that state. This ap- 
proximation reduces the dimensions of the Hamiltonian 
blocks Hab, Eq. J2J, and allows for a perturbative solu- 
tion. 

Since the only available DQD states in the strong 
bias regime are the triplet |Tq) and two singlets, \S) 
and \Dj), the Htt block of Eq. is one-dimensional, 
Hst and Hts are two-dimensional, and Hss is four- 
dimensional. For the present discussion of the strong 
bias regime, we choose the zero of the energy scale at 
ATilu -2e + U + 2V+ + U N + 2U M + U F , setting the expec- 
tation value of the energy of four singly occupied QDs 



with the DQDs in the electron singlet states to zero, 
(S, S\H\S, S) — 0. Using the expressions for the Hamil- 
tonian matrix elements given in the Appendix A, we find 
the matrices of the H a b blocks (a, b = S,T). The energy 
of the \TT) state is then 



E TT = 2(V- - V+ 



(14) 



The two-dimensional blocks Hts and Hst are related 
by the symmetry under exchange of the double quantum 
dots L <-> R and in the bases {\S,T ),\Di,T }}, and 
{|To,j§), |Tq, Dj}}, have the identical matrix form 



Hts — Hst — V. 



V2t s 
V2t s V D -e 



(15) 



where ts = —t^ + Ts is the renormalized hopping matrix 
element and Vd = U — V+ + Un — Up is the electrostatic 
energy cost of doubly occupying the pk = I state in the 
presence of the triplet DQD. With our choice of the zero 
of the energy scale, the ground state energies of Hst and 
Hts are 



Est = Ets 



V- - V + + l - (V D - e) - \ J (V D - + 8t|. 



(16) 



From the energies Est and Ets, we extract the 
isotropic exchange part of the low-energy four-spin 
Hamiltonian Eq. as 

J = Ett — Hst — Ett — Hts- (1^) 
The resulting exchange interaction strength is 



1 



1 



J = V- - V+ - - (V D - e) + -^(V D - e) 2 + 8t|. (18) 

Comparing this result with the case of an unbiased iso- 
lated double quantum dot, Eq. I|10l) . we see that the ef- 
fect of the strong bias e and the presence of another DQD 
behind the impenetrable barrier is the change of the ef- 
fective on-site repulsion to the value Vd — £ and a reduc- 
tion of the effective tunneling matrix element because 
of the large gap to the excited doubly occupied state. 
As a consequence of this gap, the isotropic exchange in 
the limit of noninteracting DQDs and weak tunneling is 
J = V- — V + + 2t 2 H /(U — V + — e), with the hopping contri- 
bution reduced to a half of the result expected from the 
standard Hubbard model in the unbiased case, 4t 2 H /UH 
@. 

The four-dimensional block Hss in the basis 
{|5,5),(|5,D/) + \D I ,S))/^/2,\D I ,D I ),(\S,D I ) - 
\D x ,S))/y/2,}]& 



Hss = 








2t s 




2t s 


V D -e + 




2X D 


2ti 











2Xc 



2X D 
2ti 
Edd 




V D 



21c 



(19) 
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where </ is the tunneling matrix element renormalized by 
the spectator DQD in the doubly occupied state, and 

E DD = 2U + W N - 2U M -U f - 2V+ - 2e, (20) 

accounts for the repulsion energy of four electrons in the 
Pk = I orbitals and the bias e, see Appendix A. Due 
to the symmetry with respect to exchange of the DQDs, 
L R, the antisymmetric state (\S,Di) — \Dj,S))/\/2 
decouples from the other, symmetric, states. 

In the limit of large and positive bias, |e — Vd\ ^> 
tg/i, Xg/D, all the tunnelling and correlated hopping 
terms in the Hamiltonian Hss can be taken to be small. 
The unperturbed Hamiltonian is then diagonal and the 
ground state energy is Edjj. This situation is relevant, 
because all the small terms are proportional to the over- 
lap S of the localized states in the quantum dots, which 
is small for weakly tunnel-coupled QDs, and we can reach 
this regime by applying external voltage to make |e — Vh\ 
large enough. 

Operating the system in the strong bias regime causes 
a qualitative change to the effective low-energy Hamilto- 
nian by turning on the entanglement generating term E e 
in Eq. ©, 

E e = Ett — 2Est + Ess- (21) 

For weak bias and in the absence of tunneling, the en- 
tanglement generating E e term is zero, as can be checked 
from the energies of the states \S, S), \Tq, S), and \T , To), 
given in Appendix A. This is not true in the case of strong 
bias, where the entangling interaction of the strength 
E e — Un~ 2Um + Up 7^ is present even if the tunneling 
terms are zero. In the strong bias regime, the conditions 
for E e — are ti = ts, Xs — X D , and E DD — 2(Vd — e). 
While the first two conditions are satisfied when there is 
no tunneling, the third is independent of the tunneling. 
It is only satisfied in the limit of long distance between 
DQDs, a-2 ^> oi, see Fig. ^ The tunneling causes a 
second-order correction to Ess, 

Ess = E DD H — h (22) 

Edd — (Vb — e) Edd 

and the corresponding correction to E e |28[ . 

We have calculated the matrix elements of the 
Coulomb interaction using the basis of single-electron 
Wannier states obtained by orthogonalizing the Fock- 
Darwin ground states centered at the quantum dots po- 
sitions, following 0. The resulting matrix elements can 
all be expressed in terms of the distances between the 
quantum dots, and the tunneling matrix element t be- 
tween QD in DQD. These results are summarized in Ap- 
pendix A. Together with Eqs. ©, dJl, and they 
provide a model of the low-energy Hamiltonian of a pair 
of qubits realized on a DDQD in the strong bias regime. 
This model can describe a two-qubit quantum gate real- 
ized by adiabatically switching the value of the control 
parameter e so that the qubit goes from the weak bias 
regime to the strong bias regime and back. 



In summary, the interaction of the DQDs causes a 
change in the parameters of the extended Hubbard model 
coupling strength, Eq. IjlUI) , so that the energies and hop- 
ping matrix elements on one DQD depend on the state 
of the other. Also, the processes in which the hopping 
of electrons on the two DQDs is correlated and medi- 
ated by the direct Coulomb interaction become possible, 
see Fig. |21 The coupling between the DQDs causes an 
effective spin interaction that deviates from the form of 
exchanged-coupled qubits, adding the entangling term E e 
to the Eq. ©. This deviation creates the entanglement 
between the two qubits. The generation of entanglement 
can be efficiently controlled by changing the bias e. 



IV. GENERAL BIAS 

The study of a double double quantum dot (DDQD) 
system in the strong bias regime presented in Sec. Ill 
allows for a simple pcrturbative solution and offers an 
insight into the mechanism of entanglement generation. 
However, it lacks sufficient predictive power for a general 
analysis of a realistic two-qubit quantum gate: When 
switching on and off the entangling interaction, a contin- 
uous voltage pulse is applied, and the system undergoes a 
smooth transition from the strong bias regime to the un- 
biased (or merely biased) regime and vice versa. During 
this transition, the system has to pass through an inter- 
mediate weak-bias regime where the perturbative expan- 
sion Eq. (|22|l breaks down. 

In this section, we calculate the full Hund-Mulliken 
(HM) Hamiltonian of the four quantum dots, including 
both \Dj) and \Do) states. This calculation allows us 
to predict the quantum gate generated by an arbitrarily 
shaped adiabatic pulse of the control parameters t and 
e. The main difference in the system's description is that 
now we take into account both doubly occupied states 
\Dj) and \Do) in each DQD. Therefore, we are working 
in the entire Hilbert space of the HM approximation, and 
the strong bias requirement is not important. Now, Htt 
is one-dimensional, Hst and Hts are three-dimensional, 
and Hss is nine-dimensional. 

Following the discussion of Sec. Ill, the effective low- 
energy spin Hamiltonian H, Eq. @, is determined by 
the energies E a b, where a,b — S, T, of the lowest energy 
states of a given spin configuration. Due to the L <-> R 
symmetry, H is the sum of the isotropic exchange terms 
and the entangling term. We proceed by calculating the 
matrix elements of the Hamiltonian as a function of the 
tunneling matrix element t and the bias e. The results 
of this calculation are given in Appendix A. Numerical 
diagonalization of the resulting Hamiltonian gives the en- 
ergies E a b, for each of the blocks H a b, where a,b = S, T. 
Finally, we extract the effective low-energy Hamiltonian 
parameters J and E e using Eq. (|17f) and Eq. I|21|l . 

The dependence of the isotropic exchange coupling on 
the bias J(e) is illustrated in Fig. EI ln the zero-tunneling 
limit, we can identify three regions of qualitatively dif- 
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FIG. 6: Isotropic exchange coupling J as a function of the 
bias e. In the regions of strong positive and negative bias, 
the exchange coupling is approximately linear J oc |e|. In the 
intermediate region, the exchange is zero in the zero tunneling 
limit and becomes nonzero as the tunneling is turned on. The 
coupling J is always positive in the absence of a magnetic field. 
The external magnetic field drives J to negative values in a 
relatively wide range of values of the tunneling matrix element 
and bias. The confinement energy of the quantum dots is 
chosen to be hu>o — 3meV, which corresponds to a quantum 
dot Bohr radius ae = 20 nm in GaAs. The distances between 
the dots are chosen to be 2ai = 1.6 ob and 2a2 = 3aB- 



FIG. 7: Entangling interaction E a as a function of bias. The 
plots correspond to different values of the tunneling matrix 
elements t within the double quantum dots in the absence 
of a magnetic field and in an external magnetic field of B = 
2 T. The t = plot indicates the regions of different lowest 
energy singlets and the positions of crossings. The strength 
of the entangling interaction E a can be changed significantly 
by tuning the bias e at a fixed tunneling matrix element t. 
Parameters used in this plot are the same as in Fig. |^1 



ferent behavior of J(e). For strong and negative bias, 
£<[/_, corresponding to the \DoDo) lowest energy 
singlet state, the isotropic exchange coupling is decreas- 
ing linearly with the bias. In the intermediate region 
U- < e < U + the exchange coupling is absent. For strong 
and positive bias [/+ < e, the exchange coupling grows 
linearly with e. The asymmetric placement of the J = 
plateau is a consequence of the different repulsion ener- 
gies of the electrons in the inner and outer QDs. As the 
tunneling is turned on, the isotropic exchange couplings 
becomes larger due to the mixing of the doubly occupied 
states in the plateau region. For zero magnetic field, the 
coupling J is positive. In a finite field there is a region 
with negative J, consistent with the analysis of [t| and 
the experimental findings of p9l |. 

A plot of the entanglement generating interaction E e 
is given in Fig. [7| The zero-tunneling value of E e shows 
a structure determined by the Coulomb energies of the 
basis states Eqs. Q-Q. In a wide plateau of small bias 
the entangling interaction vanishes, because all of the 
lowest-energy states of definite spin are products of \S) 
and \Tq). Since the direct exchange interaction V- — V+ 
is zero in the absence of tunneling, those two states are 
equal in energy. When the bias overcomes the on-site re- 
pulsion, the lowest energy states of Hss, H$t, and Hts 
change. The degenerate lowest energy states of Hss are 
either \SDj) and \DjS) in the region of large bias on the 
right of the plateau, or \SDq) and \DqS) in the region 
of smaller bias to the left of the plateau. Simultaneously, 
the analogous states with \S) replaced by \Tq) become 



the lowest energy states in Hst and Hts- In these two 
regions E e is a linear function of e, E e = Un — Uf — U — e 
on the left and E e = —Un + Uf — U + e on the right of 
the plateau. When the absolute value of the bias is even 
higher, the lowest energy state in Hss is \DjDi) for a 
very strong and positive bias and \DoDo) for a very 
strong and negative bias. These regions are character- 
ized by an e-independent E e = Un — 2Um + Uf for large 
|e|. The values U± for the bias e at which the changes in 
zero-tunneling lowest energy states occur depend on the 
geometry of the device, described by the distances 2a\ 
and 2d2 (Fig. ^) and the quantization energy Tiloq, and 
correspond to the changes in behavior of the exchange 
coupling strength J. 

The zero-tunneling case shows a desirable feature in 
that E e , the quantity that determines the entanglement 
between the qubits, can be switched on and off by tuning 
e. However, the regions of different E e can not be reached 
by adiabatic pulses in the t —> limit. Turning on the 
tunneling t between the QDs will introduce transitions 
between previously disconnected regions, and the adia- 
batic gates become possible. The simple t = picture of 
the entanglement generated by a difference in Coulomb 
energies is perturbed by the transitions. It is no longer 
possible to turn off E e throughout the plateau region by 
a change in e alone. In the plateau region, E e is generi- 
cally nonzero, but small. Therefore, in order to turn off 
the entangling interaction when t is kept constant, it is 
desirable to keep t small, and to tune e to a value where 
E e = 0. 
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V. QUANTUM GATE OPERATION 

For a quantum gate applied by the time-dependent 
Hamiltonian Eq. (0) , with the parameters t and e chang- 
ing adiabatically on the time scale set by the energy 
gap between the states within the blocks H a b, the ap- 
plied gate is determined by the splittings between the 
lowest lying states in each of the subspaces of the defi- 
nite spin. If the energies of the lowest energy states in 
singlet-singlet, singlet-triplet, triplet-singlet and triplet- 
triplet subspaces are Ess(t), Esr{t) = -Erst/), and 
Errit) respectively, the gate applied by an adiabatic 
pulse starting at the time U and finishing at tf will be 
U = dia,g(<f)ss,4 > ST,4>TS> ( t>TT), with the phases 



exp — — 



E ah (t)dt. 



(23) 



With the ability to turn the entangling interaction on 
and off and perform single-qubit gates, it is possible to 
perform a CNOT gate on a pair of qubits encoded into 
spin states of DQD. We consider a quantum gate imple- 
mented by first adiabatically turning on the entangling 
interaction for a period r on , and then again adiabatically 
switching to the Hamiltonian with the entangling interac- 
tion off for the time interval r ff ■ The lowest energy states 
in each of the SS, ST, TS, and TT subspace will acquire 
a phase dependent on the control parameters e and t and 
the pulse durations. In the on state, the Hamiltonian 
that describes the ground states in all the spin subspaces 
is, up to a constant, H on — diag(£" e , J on , J on , 2 J on ), 
where E e is the strength of the entangling interaction 
in the on regime, and J on is the corresponding exchange 
coupling. After the DDQD was in the on state for the 
time r on , the applied gate is 



7 ; • 7-011 rr 

W n = exp-i—H on . 
n 



(24) 



Similarly, during the subsequent period of duration r g 
when the entangling interaction is set to zero, the applied 
gate is 



off 



■ T ° ff IT 

exp -j— H oS , 
h 



(25) 



where H g — diag(0, J g, J a s , 2 J g) in analogy with the 
on regime. The resulting gate is 



U = UosUan = exp - 



/ <j> 

A 

A 

V 2A 



(26) 



where h\ = J on T on + J ffT g is the integrated strength of 
the exchange coupling in DQD, and %<j> = E on r on is the 
integrated strength of the entangling interaction. 

The CPHASE gate, which is equivalent to CNOT up to 
single qubit rotations, is obtained when the gate param- 
eters satisfy <j> = ran and A = nir, for an odd integer m 



and an arbitrary integer n. In order to complete a CNOT, 
we follow a pulse of on-state Hamiltonian of the duration 
Ton = rmrh/E e by a pulse of the off-state Hamiltonian 
with of the duration r Q g = Ti{wk — J nT n/^)/JoS- The 
resulting gate is diag(— 1, — 1, — 1, 1 ) = -CPHASE, for 
odd n and diag(— 1, 1,1,1), which is equal to CPHASE 
with the X gate applied to both qubits before and after 
hi. For any integer n, 



CPHASE - (£®0«(£<8>0. 



(27) 



where £ = exp(i7r(l + (— l) n )a x /A). In order to complete 
the CNOT, we apply the one-qubit Haddamard gates 
H = (X + Z)/y/2 to the target qubit both before and 
after the entangling gate IA. The entire construction can 
be represented as 

CNOT = (1 <g)F)(£®0^(£ ® 0(1® H). (28) 

Note that the CNOT construction necessarily involves 
the single qubit rotations about pseudospin axes differ- 
ent from z. Such operations can be performed using the 
asymmetric bias within a DQD that encodes the qubit in 
an inhomogeneous external magnetic field |27l |. The en- 
tangling part of a CNOT gate can be performed by puls- 
ing the bias e only, and keeping the tunneling t constant. 
Therefore, control over the bias e and the availability of 
an inhomogeneous magnetic field are sufficient for the 
universal quantum computing with two-spin qubits. 



VI. CONCLUSION 

We have analyzed two-qubit gates in a pair of qubits, 
each encoded into singlet and triplet states of a DQD, 
and coupled by Coulomb repulsion. A two-qubit CNOT 
gate, which together with the single qubit rotations forms 
a universal set of quantum gates, can be performed by 
tuning the bias of the inner dots with respect to the outer 
ones. We identify the entangling interaction strength E c 
as a quantity that has to be controlled in order to imple- 
ment a CNOT with the aid of single qubit rotations. 

The dependence of E e on the externally controllable 
bias e and the tunneling matrix clement t shows that it 
can in principle be turned on and off by changing e alone, 
if sufficiently low values of t are available. 

The largest change in E c comes from tuning of the 
system through the resonance between singly occupied 
state and doubly occupied state on a DQD. At the side 
of the resonance with a singly occupied ground state, 
and far from the resonance, the entangling interaction 
E e is caused by inter-DQD correlation and is small. On 
the other side of the resonance, with a doubly occupied 
DQD ground state, the entangling interaction is caused 
by the direct Coulomb repulsion and it is much stronger. 
Two-qubit gates necessary for a universal set of gates can 
be performed by switching between the strong and weak 
entanglement generation regimes using voltage pulses. 
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APPENDIX A: HUND-MULLIKEN 16 x 16 
HAMILTONIAN 

The full Hund-Mulliken Hamiltonian is block-diagonal 
due to the symmetry of the interactions with respect to 
arbitrary rotations in spin space. In reality, this sym- 
metry is broken by the weak spin-orbit coupling inter- 
action that we have neglected. The blocks are the one- 
dimensional Htt, the two three-dimensional Hts and 
Hst, and the nine-dimensional Hss, where T stands for 
a triplet and S for a singlet state on a DQD. In this Ap- 
pendix, we present the matrices of these blocks as func- 
tions of the system geometry and the control parameters. 



Hss = 



We do not antisymmetrize with respect to the permuta- 
tions of electrons that belong to different quantum dots 
and have non-overlapping orbital wave functions. The 
matrix elements of the Hamiltonian that describe the 
Coulomb interaction within a DQD (intra-DQD terms) 
U, t, X, V + and V- were analyzed in The inter-DQD 
elements depend on the following matrix elements of the 
Coulomb interaction between the product states of the 
\qps) electrons localized in the qubit q and the quantum 
dot p, and having a spin s, 

X s = (LIs,RIs'\V c \LOs,ROs'), (A4) 

X D = {LIs,ROs'\V c \LOs,RIs'}, (A5) 

T = {LOs,ROs'\V c \LIs,ROs'}, (A6) 

T 7 = {LIs,RIs'\V c \LIs,ROs'). (A7) 

In zero magnetic field, we find that Xs = Xd- 
The off-diagonal elements are determined by 

t s = T +Tj-t„, (A8) 
ti = 2T/-t H , (A9) 
t = 2T -t H , (A10) 



There is only one TT state and its energy is 

H TT = E TT = 2VL +U N + 2U M + U f - 2e. (Al) 

The three-dimensional blocks Hts and Hst are related 
by the symmetry operation of exchanging the DQD and 
if we choose the basis {\T , S), \T , Dj), \T , D )} for the 
TS and {IS, T ), \Dj,T ), \D o ,T )} for the ST subspace, 
they can both be represented by the matrix 

/ C TS V2t s V2t s \ 
Hts = H ST = y/2t s C T i X . (A2) 
V V2t s X C TO J 

The nine-dimensional block of singlet states, in the direct 
product basis composed out of the two-electron states 
Eq. is 



(A3) 



I 

and the diagonal elements are given by 



Ctt 


= 2V-- 


f U N + 2U M + U F 


-2e, 


(All) 


Cts 


= V+ + 


V-+U N + 


2U M 4- 


U F - 2e, 


(A12) 


Cti 


= V- + 


U + 2U N + 


2U M - 


3e, 


(A13) 


Cto 


= V- + 


U + 2U M + 


2U F - 


e, 


(A14) 


Css 


= u N ^ 


2U M + U F 


+ 2V+ 


-2c, 


(A15) 


Csi 


= 2U M 


+ 2U F + U 


+ V+- 


3e, 


(A16) 


Cso 


= 2U M 


+ 2U F + U 


+ V+- 


c, 


(A17) 


Cn 


= w N 


+ 2U~4e, 






(A18) 


Cio 


= W M 


+ 2U-2e, 






(A19) 


Coo 


= AU F - 


f 2U, 






(A20) 



where the symmetry with respect to exchange of the 
DQDs leads to C A b = C B a where A,B e {T, S, I, O}. 

To represent the matrix elements in terms of the sys- 
tem parameters, the single QD quantization energy Hujo, 
tunneling matrix element within an isolated DQD t, the 
bias e and the interdot distances a\ and <22, we have 
to adopt a model for the binding potential of a DQD 
and the orbitals of Hund-Mulliken approximation. We 
assume that the QD orbitals are Wannier functions ob- 



J 



/ Css 


V2t s 


V2t s 


V2t s 


2X D 


2X S 


V2t s 


2X S 


2X D \ 


V2t s 


Csi 


X 


2X S 


V2ti 





2X D 


V2~h 





y/2t s 


X 


Cso 


2X D 





V2t 


2X S 





V2t 


V2t s 


2X S 


2X D 


Cis 


y/2ti 


V2tj 


X 








2X D 


V2tj 





V2t! 


Cn 


X 





X 





2X S 





V2t 


V2tj 


X 


Cio 








X 


V2t s 


2X D 


2X S 


X 








Cos 


V2t 


V2t 


2X S 


y/2ti 








X 





V2t 


Coi 


X 


\ 2X D 





V2t Q 








X 


V2t Q 


X 


Coo / 



11 



tained by orthogonalization of the Fock-Darwin ground 
states centered at the positions of the QDs within a DQD, 
(a 2 ± a±, 0) and (— a 2 ± ax, 0). The Wannier orbitals are 
of the generic form 

\W qJ ) = N (\(j) qt i) — g\cp q ,o)) j (A21) 
|W g ,o) =^(-g|0 9 ,/> + l^o», (A22) 

where \4> q j(o)) is the Fock-Darwin ground state on the 
dot belonging to the qubit q = L, R and the inner(J) or 

I 

U N = cN 4 

AgS 
U F = cN 4 

AgS 

U M = cN 4 
2gS\ 
T = cN 4 

(9 + 
T) = cN 4 

(9 + 

X s = cN 4 

2S (< 
X D = cN 4 
25 (f 

in terms of the overlaps of the harmonic oscillator wave 
functions S, the mixing factor g, and the function 

f(d, 1) = V^exp ( — ot{d, I)) I Q (a(d, I)) , (A30) 

where a(d, I) = bd 2 — (b— l/b)l 2 . We use the contraction 
factor b = lu/loq to measure the magnetic field strength. 
The o veral l strength of the Coulomb interaction is set by 
c = \/ 7r/2e 2 / nhajQCiB: where e is the electron charge, k 



outer (O) QD, Eq. (JT2J. The Wannier orbitals are de- 
termined by the overlap of these wave functions, S = 
(0gj|<A?,o) = exp (— d\ (2b — 1/b)), through the mix- 
ing g = (1 — yl — S 2 ) j S and normalization constant 
N = l/^l-2gS + g 2 . 

The Coulomb interaction matrix elements for the DQD 
centered at ±02 = i^ae and QDs within a DQD dis- 
placed by ±<zi = ±diaB from the center of the DQD are 
then expressed as 



(A23) 



(A24) 



(A25) 



(A26) 



(A27) 



(A28) 



(A29) 



is the dielectric constant, and hu>o is the single isolated 
QD quantization energy 

To model the dependence of the matrix elements on 
externally controllable tunneling matrix element t, we 
use the connection between the tunneling and the overlap 
S = S(t) that holds for the quartic double well, Eq. i|13|) 
and assume that it holds throughout the gate operation. 



/(da - dx, 0) + 2g 2 (l + S 2 ) /(d 2 , 0) + g 4 f(d 1 + d 2 , 0) + 2S 2 g 2 f(d 2 ,d 1 )~ 
/(da + dt, 0) + 2 3 2 (1 + S 2 ) /(da, 0) + 5 4 /(d 2 - di, 0) + 2S 2 g 2 f(d 2 ,dt)- 

( /(4+ ** )+ ^ ( * _**));, 

(1 + , 9 4 )/(d 2 , 0) + g 2 (/(dx + d 2 , 0) + /(di - d 2 , 0) + 2S 2 (/(d 2 , 0) + /(da, dr))) - 

S ((1 + 3.9 2 ) /(d 2 + I , I ) + ( 5 4 + 3. 9 2 ) /(d 2 - y , y )) - 
5 3 ) ((1 + S 2 ) /(d 2 , 0) + S 2 .f(d 2 ,d 1 )) - 3 /(d 2 + d l5 0) - 5 3 /(d 2 - dx, 0) , 
5 ((1 + 3g 2 ) f{d 2 y , y ) + ( 5 4 + 3. 9 2 ) /(d 2 + y , y )) - 
5 3 ) ((1 + S 2 ) /(da, 0) + S 2 f(d 2 ,dx)) - gf(d 2 - dx, 0) - 5 3 /(d 2 + d lf 0) , 
(S 2 + 2. 9 2 + , 9 4 5 2 ) /(d 2 , 0) + 3 2 (/(da + d 2 , 0) + /(d a - d 2 , 0) + 2S 2 /(d 2 , di)) - 

, + 9 », ( /((i3 + **, _**,)■, 

S 2 (1 + ff 4 ) /(da, di) + ff 2 (/(di + da, 0) + /(d 2 - d 1; 0) + 2 (l + S 2 ) /(da, 0)) - 

s\ f , f , di di . , di di A 
? + 3 3 ) ( v /(d 2 + y , y ) + /(rf 2 - y , y )J , 
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